Structural and dynamical properties of a gMasz-one-dimensional classical binary system 
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The ground state configurations and the "normal" mode spectra of a quasi-one-dimensional (QID) 
binary system of charged particles interacting through a screened Coulomb potential are presented. 
The minimum energy configurations were obtained analytically and independently through molec- 
ular dynamic simulations. A rich variety of ordered structures were found as a function of the 
screening parameter, the particle density, and the ratio between the charges of the distinct types of 
particles. Continuous and discontinuous structural transitions, as well as an unexpected symmetry 
breaking in the charge distribution are observed when the density of the system is changed. For 
near equal charges we found a disordered phase where a mixing of the two types of particles occurs. 
The phonon dispersion curves were calculated within the harmonic approximation for the one- and 
two-chain structures. 

PACS numbers: 64.60.Cn, 82.70.Dd, 63.20.-e, 83.10.-y 
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I. INTRODUCTION 

Advances in experimental and numerical techniques al- 
lowed a considerable improvement of the possibilities to 
create and study new crystalline structures which will 
open a vast field of scientific and technological appli- 
cations. Recently, a variety of novel three-dimensional 
(3D) binary super-lattices obtained through the combina- 
tion of semiconducting, metallic, and magnetic colloidal 
nano-particles was observed experimentallyi'=. The self- 
organization in a two-dimensional (2D) binary colloidal 
system resulted in interesting mixed configurations with 
partial clustering"^ . The use of colloidal particles is tech- 
nically interesting since it allows a real time and spatial 
direct observation of their position which is a great ad- 
vantage as compared to atoms or molecules, as shown 
recently in an interesting experimental study concerning 
defect induced melting^. 

Systems of particles moving in space of a reduced di- 
mensionality or submitted to an external confinement ex- 
hibit different behavior from their free-of-border coun- 
terparts. It has been shown previously that the pres- 
ence of constraints generates new features which are in 
some cases unexpected, as e.g. the curious ground state 
configuration found in a 3D finite size binary Coulomb 
cluster- , where a mixed arrangement of the two types of 
particles was observed in the core, while the edge of the 
system was characterized by shell structures of separated 
species. Similarly, in 2D systems we can cite a reentrant 
melting phenomenon in a 2D circular confined cluster of 
repulsive dipoles^ii, and interesting asymmetrical ground 
state configurations observed in symmetrically confined 
Coulomb clusters^i^. 

Under specific conditions of density and temperature 
an electron gas crystallizes in an ordered structure form- 
ing the so called Wigner crystal (WC)^°. In general, 
a Wigner crystal structure can also refer to an ordered 
phase occurring in other systems of interacting particles. 



Several non-electronic systems show such a Wigner or- 
dered phase, which has been experimentally observed 
in colloidal systemaii"— , dusty plasmas^^"^^, and ion 
trapsiliiS. 

Depending on the dimensionality and extent of the sys- 
tem the ordered phase can have different symmetry. As 
shown theoretically, the minimum energy arrangement 
of the 2D infinite electron gas is the triangular latticei^, 
while the circular confined analog cluster was observed 
experimentally^ and numerically^!, and has a mixed 
ground state configuration with a triangular structure 
in the core and a circular structure at the edge of the 
system as a consequence of the symmetry of the external 
confinement^SiSi. 

When the 2D infinite system is subjected to an extra 
confinement (e.g. parabolic) in one direction the system 
can be called guasi-one-dimensional (QID). Such a QID 
system of charged particles interacting through a repul- 
sive potential self-organizes in a chain-like structure were 
recently experimentally studied in dusty plasma o^^d^ and 
paramagnetic colloids^^i^, and thr oug h analytical and 
numerical calculations^^i^. In Ref. [2J, the authors pre- 
sented a systematic and interesting study of the struc- 
tural, dynamical, and thermal properties as a function of 
the linear density of the system. It is claimed that the ID 
Coulomb crystal confined in a storage ring may be used 
in atomic clocks or even for quantum computation^^"^^. 
An interesting application for the columnar regular struc- 
ture of interacting particles was shown in Ref. l29l where 
a regular linear array of superparamagnetic colloidal par- 
ticles confined in a thin gap was used for the separation 
of DNA. 

Besides the single-component system of interacting 
particles, which has been widely studied in the last years, 
Wigner crystallization may also be observed in multi- 
component systems of interacting particles^fl. The sim- 
plest multi-component case is the binary mixture of two 
types of particles which exhibits a richer set of physi- 
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cal properties when compared to the one-component sys- 
tem. For example, a large number of different equilib- 
rium configurations which are intrinsically dependent on 
the relative fraction of the different types of particles was 
recently found in a 2D system of dipoles'^^. An interest- 
ing analysis of the pair correlation function decay as a 
function of the packing fraction in a binary system of 
hard-spheres was presented in Ref. ,32> The presence 
of particles with distinct physical properties (e.g. size, 
charge, mass) introduces a competition between different 
scales, which is the reason for the richer phenomenology 
in such systems. 

In the present paper, we extend the work of Ref. [U 
to a binary two-dimensional system of charged particles 
interacting through a screened Coulomb potential and 
confined in a parabolic channel. The particles are al- 
lowed to move (without friction) in the x-y plane, but due 
to the external parabolic confinement in the y-direction, 
the system acquires a QID character, i. e. the parti- 
cles are free to move in the x-direction and have a con- 
strained motion in the y-direction. Due to the possibility 
to tune the range of the screened interaction potential in 
the present model, the binary system of hard-spheres'^^ 
becomes a particular case of our system, which occurs 
when the screening length of the interaction potential is 
very small. In spite of the specific interaction potential 
considered in the present paper, the main qualitative fea- 
tures of the obtained results are expected not to depend 
on the functional form of the inter-particle interaction. 

We present a systematic study of the structure of the 
ground state as a function of the linear density, the ra- 
tio between the charges of the distinct types of particles, 
which characterizes the binary system, and the param- 
eter which determines the range of the repulsive inter- 
action potential between the particles. The range of the 
interaction between the charged particles can be varied in 
experimental systems of colloidal particles by changing, 
for instance, the salt concentration in the solvent^^. In 
the present analysis we limit ourselves to an equal den- 
sity of the two types of particles which are assumed to 
have the same mass in order to reduce the number of 
parameters. The normal modes, i.e. the phonons, of the 
present system are also studied. 

The paper is organized as follow. In Sec. II, we de- 
scribe the model and the procedure used to find the min- 
imum energy configurations and the frequencies of the 
normal mode spectra. In Sec. Ill we study two different 
zero temperature phase diagrams. In the first, the ground 
state configurations are presented for distinct screening 
parameter of the interaction potential and densities. In 
the second, the dependence of the crystal structures are 
shown for different values of the ratio between charges 
and density. The analysis of the phonon dispersion curves 
is presented in Sec. IV. Our conclusions are given in Sec. 
V. 



II. NUMERICAL APPROACH 

We studied a two-dimensional binary cluster consisting 
of a large and equal number of particles with distinct 
charges qa and qb, which are allowed to move in the x-y 
plane. The charged particles interact through a repulsive 
Debye-Hiickel (or Yukawa) potential exp{—r/X)/r, are 
free to move in the a;-direction, and are confined by a one- 
dimensional parabolic potential which limits the motion 
of the particles in the y-direction. The potential energy 
of the system is given by 
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where e is the dielectric constant of the medium the par- 
ticles are moving in, A is the Debye length, and = |ri| 
is the distance of the z"* particle from the center of the 
confinement potential. In order to reveal the important 
parameters of the system, it is convenient to write the en- 
ergy and the distances in units of Eq = {m(jjQq'^/2e^y^'^ 
and To = (2g^/mea;g)^/'^, respectively, and to define the 
quantity a = qt/qa (with qa = e is the reference charge 
taken to be equal to the elementary charge), and the 
screening parameter k = ro/A. In so doing, the expres- 
sion of the potential energy is reduced to 
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and the state of the system is determined now by the 
ratio between charges a, the number of particles (which 
can be associated to the density), and the dimensionless 
screening length k. The temperature is expressed in units 
of To = Eo/ks, where ks is the Boltzmann constant. 

The minimum energy configurations were found on the 
one hand by analytical calculations, and on the other 
hand also through molecular dynamics simulations. In 
the numerical simulations we typically considered 300 
particles, together with periodic boundary conditions in 
the unconfined direction in order to mimic an infinite 
system. 

The present model system does not address effects due 
to frictional dissipation. In spite of the primary impor- 
tance of friction to the motion of the particles in real 
systems, the ground state configurations are not affected 
by it. 
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FIG. 1: (Color online) Ground state configurations. Tlie y = axis is shown by the thin horizontal dotted line and the two 
different types of particles in the binary system are indicated by different symbols. The red open circles are particles with 
charge qb/qa = a = 2, and the black solid spheres are particles with qa/e = 1. The linear density of each configuration is also 
shown. 



III. GROUND STATE CRYSTAL STRUCTURES 

In view of the large parameter space, and in order to 
show the increased richness of a binary system we discuss 
first the case a = 2. We only considered equal mass 
particles. If masses are taken different the two types of 
particles will feel a different confinement potential which 
will result in a quantitative modification of the phase 
diagram^. 

In this section we will present the structural results for 
our binary QfD system. The results were obtained both 
analytically and numerically. In the former, we calculate 
the energy per particle for configurations with different 
number of chains, different arrangement of the particles 



within the chains, and minimize such expressions with 
respect to the distance between chains. To study the 
dynamical properties of the system we use the harmonic 
approximation in order to obtain the normal modes of 
the ground state configuration. This is shown in section 

llYl 

Before starting to present our results it is convenient to 
introduce here a dimensionless linear density (n) , defined 
as the ratio between the total number of particles in the 
unit cell of the crystalline arrangement and the length of 
the unit cell in the unconfined direction. This is similar 
to the linear density defined previously in Ref. [24). 

First, we present the minimum energy configurations 
and structural transitions as a function of the range of the 
interaction potential and the linear density of the system. 
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FIG. 2: (Color online) (a) The energy per particle, and (b) 
the lateral position of the chains as a function of the density 
for K = 1 and a = 2. The black solid circles are the particles 
with qa/e. = 1 and the red open circles with q^/e = 2. (c) The 
first derivative of the energy with respect to the density. 



More specifically, we study how the ground state configu- 
rations depend on the the screening parameter (k). Next, 
we study the system as a function of the ratio between the 
charges of the two types of particles (a). In the present 
study we limit ourselves to an equal fraction of particles 
with distinct charges. 



A. Dependence on k 

In general, the different charged particles organize 
themselves in chains. The number of chains, and how 
particles are placed in them, depends on n , the screening 
parameter (k), and the ratio between the charges (a). In 
this sub-section, we show how the crystalline structure of 
the QID system depends on k and n. The ratio between 
charges will be kept fixed a — 2.0. 

Lets first consider the case with k = 1. When n in- 
creases the following sequence for the first five ground 
state configurations is observed: one chain — two chains 
— )■ four chains (case 1) — four chains (case 2) — )■ seven 



chains (see Fig. [T]). The energy per particle as a function 
of n is shown in Fig. [H^a). As can be seen in Fig. [TJ in the 
one-chain regime the particles with distinct charges are 
arranged alternately, and equally spaced on the a:;-axis, 
where the confinement potential is zero. In this case, the 
unit cell consists of two particles, one of each type (Fig. 
[T]). The linear density is n = rg/a, and the x coordinate 
of the particles with distinct charges are Xi = 2ia and 
Xi — {2i — l)a, with i — ±1, ±2, ±3, .... The energy per 
particle for the one-chain regime is 
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were the first term accounts for the interaction between 
particles of the same type, while the last one represents 
the interaction between particles with different charges. 

When n ~ 0.74155 the two different particles start to 
move off chain, alternating in the positive and negative 
y-direction, leading to a new ground state configuration 
consisting of two chains. We obtained the following ex- 
pression for the energy per particle. 
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(4) 



where c — d/a is the dimensionless distance between the 
chains (Fig. [1]). The two-chain regime is characterized 
by a remarkable asymmetry in the charge distribution. 
It is interesting to comment that a similar behavior was 
experimentally observed previously in a two-component 
Coulomb crystals in a linear (quasiSD) Paul traps and 
in a binary system of charged dry grains^i^. The charge 
segregation characterizing the two-chain regime is a con- 
sequence of a spontaneous symmetry breaking which oc- 
curs through a continuous structural transition. This is 
made clear in Fig. [2l^b), where the lateral position of 
the chains is plotted as a function of the density. We 
found that with the exception of the transition from the 
one-chain to the two-chain regime, all other structural 
transitions are characterized by a discontinuity in the 
lateral position of the chains which is associated with a 
first order "structural" transition. In this case, the first 
derivative of the energy with respect to the density is 
discontinuous at the transition point [Fig. [2][c)]. For the 
continuous transition (second order structural changes) 
a discontinuity is observed only in the second derivative 
of the energy with respect to the density. 

A further increase in the density brings the system to 
the four-chain regime. In this case, the different types 
of particles are again symmetrically distributed with re- 
spect to the X-axis, but still segregated. The four-chain 
configuration is observed in two different minimum en- 
ergy configurations, defined here as case 1 and case 2. In 
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FIG. 3: The zero-temperature K — n phase diagram for a = 2. 



the four-chain (case 1) regime the rows with the same 
type of charges are displaced with respect to each other 
by a distance a/4 along the unconfined direction (x-eods), 
while neighbor rows with distinct types of particles are 
displaced by a distance a/2 along the a;-axis (Fig. [1]). 
Also, the distance between the internal chains (consist- 
ing of particles with the same charge) is larger than the 
distance between the internal chains and the external 
ones [Fig. [Hb)]. This is interesting because the inter- 
action between chains with distinct charges is intuitively 
expected to be larger than the interaction between chains 
with the same lower charge. In the four-chain (case 2) 
regime, chains with the same charge are displaced by a/ 2 
with respect to each other along the unconfined direction. 
The distance between chains has an opposite behavior to 
the one of case 1, i.e. the distance between chains with 
distinct charges is larger than the distance between the 
internal rows [Fig. [1Kb)]. 

The specific expressions for the energy per particle for 
the four-chain configurations, as well as for other arrange- 
ments presented in Fig. [I] can be obtained directly from 
Eq. (dJ. The three-chain configuration shown in Fig. 
[1] does not become a minimum energy arrangement for 
a = 2. However, as will be shown in the next section, 
for some particular values of a such a configuration was 
found as the ground state. 

As can be observed in Fig. when the density is in- 
creased the system tends to crystallize in structures con- 
sisting of a large number of chains. In such cases it is not 
trivial to predict the minimum energy arrangement of 
the particles. To guide our analytical analysis we resort 
to numerical molecular dynamic simulations. Due to the 
larger number of particles considered in the simulations, 
there are many stable configurations with energy only 
slightly different from the minimal one. In fact, it is ex- 
tremely difficult to reach the ground state configuration, 
even in our numerical simulations. The final configura- 
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FIG. 4: (Color online) (a) The zero-temperature a — n phase 
diagram for k = 1. (b) Zoom of the small a and intermediate 
density regime. Solid (dotted) lines are first (second) order 
structural transitions. The symbols are the calculated points. 
The red dashed region corresponds to the values of a and n 
for which the two types of particles can be found mixed in the 
same chain. 



tion typically observed in the simulation was a mixture 
of several structures. We used this output as a hint to 
search for the minimal energy arrangement. In Fig. [1] 
only the stable configurations that have minimal energy 
are presented, but many others were considered. 

The minimum energy structures presented so far were 
also found when we minimized the analytical expressions 
for the energy per particle with respect to the distance 
between chains for different values of k and n (a = 2). 
The result is presented as a zero-temperature n — n phase 
diagram in Fig. [31 As a general feature, the interval of 
density in which a particular phase is observed increases 
with increasing n. Note that the larger the value of k 
the shorter the range of inter-particle interaction poten- 
tial. This means that the critical distance between adja- 
cent particles necessary to produce an interaction strong 
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FIG. 5: (Color online) Example of the mixed configuration 
obtained from molecular dynamic simulations. Black solid 
circles represent the reference charge Qa/e = 1, while red open 
symbols are the other particles with charge qb/e = a. 



tern presents the following sequence of five different min- 
imum energy structures as the density is increased: one 
chain — )• two chains — >■ four chains — )• three chains — > 
four chains five chains. These transitions occur re- 
spectively at densities n=0.8946, 2.0312, 2.1389, 3.3222, 
4.7054. Notice that such a sequence is different from the 
one found in the case a 1 [Fig. IUJa)] . It should be em- 
phasized that the three-chain structure observed when 
a = 1 is distinct from the one presented in Fig. [TJ In the 
three-chain configuration of the equal charged particles 
case the rows present the same density and are displaced 
with respect to each other by a distance a/2 along the x- 
axis. This is not the case for the three-chain arrangement 
commented before when 0.136 < a < 0.237 for interme- 
diate values of density [see Fig. [Tj. In the next section, 
we directed our attention to the structure of the present 
binary system when the distinct charges are almost equal 
(««1). 



enough to change the phase structure will be smaller, and 
consequently, the critical density will be larger. 



B. Dependence on a 

Now we will study how the crystal structure of our 
system depends on the density (n) when different values 
of the ratio between charges (a) are considered. In this 
section, we fixed the screening parameter to k = 1 which 
is a typical value in colloidal systems and dusty plasmas. 

For different values of a the configurations are the same 
as the ones presented in Fig. [H but which type of particle 
is located in the internal and in the external chains de- 
pends whether a is < 1 or > 1. In general, particles with 
smaller charge are located in the more internal chains. 

The zero temperature a — n phase diagram is presented 
in Fig. mja). For a > 1, the same sequence of structures 
previously shown in Fig. [3] are found. The behavior is 
very different for a < 1. In this case, new structures 
appear or disappear depending on the values of a and 
n. For example, an arrangement with three chains (see 
Fig. [1]) is observed for a < 0.237, as seen in Fig. (Hb). 
In the interval 0.136 ^ a < 0.237 the three-chain regime 
appear as the ground state configuration for intermediate 
density. The two-chain — ^ three-chain transition is a first 
order structural transition. 

For a < 0.136 the two-chain regime is no longer ob- 
served as the ground state for any value of n. Instead, 
the system changes continuously (second order structural 
change) from the one- to the three-chain configuration 
shown in Fig. [1] as the density is increased [Fig. |H[b)] 
through a zig-zag transition. The transition three-chain 
— > four-chain (case 1) is a first order transition for all 
values of a. 

For a — 1 the particles are identical, and the present 
model is identical to the one studied previously in Ref. 
[2^ . We found, in agreement with Ref. [l^, that the sys- 



C. Region around a — 1: disordered phase 

Notice that in the density region I ^ n < 2 the 
four-chain (case 1) phase is strongly suppressed in fa- 
vor of the two-chain phase, which for a = 1 extends 
now up to n = 2.0312. Beyond this density, there is 
a small density interval where the four-chain (case 2) 
regime is found which for a = 1 is identical to the pre- 
viously found four-chain phased. In the density interval 
2.0312 < n < 3.3222 we found a small region [red shaded 
area in FiglDJa)] where the two types of particles are al- 
most randomly distributed over three-chains which for 
a ^ I were found to be no longer perfectly straight (see 

Fig. El). 

This disordered region, where the chain-like structures 
are the same as the ones observed in the case a = 1, are 
characterized by the presence of particles with distinct 
charges at the same row. Note that the red dashed region 
of the a~-n phase diagram is inside the region associated 
with the four-chain (case 2) regime [Fig. [DJa)]. For n > 
3.3221 the four-chain (case 2) regime [Fig. [1] is observed 
for any value of a. 

The mixed or disordered configurations with three 
chains were obtained through molecular dynamic sim- 
ulations. Only in the mixed region of Fig. [4ja), the en- 
ergy of such configurations were observed to be smaller 
than those obtained analytically from the proposed reg- 
ular structures (Fig. [T]). In Fig. [51 the percentage of 
the difference in energy between the four-chain (case 2) 
regime and the mixed configurations inside the mixed re- 
gion (dashed) of Fig. Ufa) is presented as a function of 
the density. The largest difference in energy was observed 
for the case a — 1 where it was less than 1 %. 

Note that asymmetric mixed configurations were found 
recently in a different system of binary charged parti- 
cles which were confined by a circular parabolic® or hard 
wall^i^ potential. 
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FIG. 6: (Color online) The relative difference between the 
energy of the four-chain (case 2) regime, obtained analyti- 
cally, and the energy of the mixed configuration, obtained 
from molecular dynamic simulations, for values of a in the 
mixed region presented in Fig. HJa). 
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FIG. 7: (Color online) The motion of the particles for each 
type of normal modes for the one- and two-chain regimes. 



IV. PHONON SPECTRUM 



We analyze how the normal mode spectrum of the 
present binary system (in the absence of frictional 
dissipatiortSS-) behaves as a function of a, rt, and k. Tak- 
ing into account the translational invariance of the sys- 
tem in the unconfined direction (x-axis), we calculated 
the normal modes within the harmonic approximation"^^. 
Since we are strictly dealing with a two-dimensional sys- 
tem, the number of degrees of freedom per unit cell is 
twice the number of particles in the unit cell (the unit 
cell are indicated by dotted boxes in Fig. [T|). Therefore, 
if I is the number of particles per unit cell there will be 11 
branches for the phonon dispersion curves. Half of these 
branches correspond to oscillations along the chain (x- 
axis - longitudinal modes) , while the others are associated 
with vibrations along the confinement direction (y-axis - 
transverse modes). If particles in the unit cell present an 
in-phase vibration, the mode is also dominantly acous- 
tical, while the opposite out-of-phase oscillation deter- 
mines the optical mode. In general, a normal mode can 
be classified in one of the following classes: longitudinal 
optical (LO), longitudinal acoustical (LA), transverse op- 
tical (TO), or transverse acoustical (TA). The motion of 
the particles for each type of normal mode for the one- 
and two-chain regimes is illustrated in Fig. [T] 

In the harmonic approximation the normal modes are 
obtained by solving the system of equations 



])Qv,i 



0, 



(5) 



where Qy,j is the displacement of particle j from its equi- 
librium position in the v direction, \i and v refer to the 
spatial coordinates x and y, ^fLv,ij is the unit matrix and 
^u.v,ij is the dynamical matrix, defined by 



(6) 



where u is an integer assigned to each unit cell. The force 
constants are given by 



, / ^ exp(-K^ a; - x'Y + {y - y'Y) 
^^.A'u) = d,,d, -======= , (7) 

^{x - x'Y + {y- y'Y 

with {x — x') € [au, a(u — 1/2)], and {y~y') = interchain 
distance with (x, y) and (x', y') the equilibrium positions 
of the particles in the unit cell, and 



<)fi,u{u 0) = - ^ (l>p,,i^{u). 

u#0 



(8) 



The phonon frequency is given in units of loq/^/I. As 
an example, the complete dynamical matrix for the one- 
chain regime is given in Appendix |X] as an example. 

The frequencies for the one-chain configuration are 
given by 



{Ai+A:i)±d{Ai-A:iY+AAi 



for the longitudinal modes, and by 



(9) 



ujt = \/ 1 + {A2 + Ai) ± ^{A2 - Aif -f 4Ai (10) 

for the transverse modes. The parameters Ai,A2, ■■■,Aq 
are given in Appendix|Xl The wave number k for the one- 
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FIG, 8: The phonon spectrum for the one chain configuration 
with K = 1 and a = 2 for different values of the density. 
The transverse acoustical (TA), the transverse optical (TO), 
the longitudinal optical (LO), and the longitudinal acoustical 
(LA) phonon dispersion curves are identified, respectively, by 
the solid, dashed, dotted and dot-dashed curves. 



FIG. 9: The normal modes of the system at the transition (a), 
(b), (c) from the one-chain to the two-chain regime, and (d) 
from the one-chain to the three-chain structure. The TA, TO, 
LO, and LA relation dispersion curves are given, respectively, 
by the solid, dashed, dotted, and dot-dashed curves. 



and the two-chains regimes is in units of -k /2a, where 2a 
is the length of the unit cell in the x-direction. 

In Fig, [5] the phonon spectrum of the system in the 
one-chain regime is presented for different values of the 
density n, and fixed values of k = 1 and a — 2. As can 
be observed, there is a clear dependence of the dispersion 
curves on n. For small values of the density, e.g, n = 0,2, 
the transverse modes (optical and acoustical) are almost 
the same [Fig. EJa)]. As the density increases the fre- 
quency of the TO mode decreases. An interesting feature 
to keep in mind is that specially for the TO mode, of the 
one-chain regime, there is a resulting force in the direc- 
tion of the confinement potential (y-direction) as a conse- 
quence of the repulsive interaction between particles with 
distinct charges. For a sufficiently small density, the dis- 
tance between adjacent particles becomes large enough 
to make the interaction between particles irrelevant to 
affect the oscillation of such particles around the equi- 
librium positions. In this case, the repulsive force in the 
y-direction is very small and the TO mode will be mainly 
determined by the confinement potential, A similar ex- 
planation can be given for the behavior of the TA mode, 
but in this case the repulsive force in the y-direction is 
smaller than that for the TO mode. Therefore, for a 
sufficiently small density the transverse modes will be 
determined predominantly by the restoring force due to 
the confinement potential, which is the reason why such 
modes become almost indistinguishable for small densi- 
ties. 

When the density is increased the distance between 
adjacent particles becomes smaller, the repulsive force 
between them increases and acts as a retarding force. 



As a consequence, the frequency of the TO oscillations 
decreases. 

The LO mode has the opposite behavior as compared 
to the TO mode, i.e. there is a hardening of such mode 
when the density increases. This is a consequence of the 
larger repulsion due to the closer proximity between par- 
ticles, A hardening with increasing density is also ob- 
served for the LA mode. 

For the LA mode, a linear dispersion curve as fc — > 
is found. This means that in this limit the longitudinal 
wave propagates with constant velocity along the chain, 
and such a velocity increases with increasing density (Fig, 
©, 

As can be observed in Fig, [SI specially in the limit 
/c — > 0, the TO branch decreases with increasing density. 
For n 0.74155 the system with hi — 1 and a = 2 
undergoes a second order structural phase transition to 
the two-chain regime. From Fig, [HJa) we notice that the 
TO mode is zero for fc = and becomes linear in the 
limit /c — >■ 0, similar to the LA mode. Such a behavior 
characterizes a continuous structural transition. 

In Fig, [9] the dispersion curves for systems with k = 1 
and distinct values of a are presented. In each case, we 
consider the critical density at which the system transits 
from the one-chain to the two-chain [Figs. IHta), (HUb), 
or to the three-chain regime [Figs. [9jd)]. As can 
be observed in Figs, [HIa) , Efb) , HJc) for a = 2, a = 0,3, 
a = 0.2, respectively, the TO mode is zero for k ~ and 
is linear for small fc-values. 

In Fig, injd), a = 0,1, the dispersion curve for the TO 
mode is qualitatively different from the previous ones. 
The softening of the TO mode is still observed, but now 
at the edge of the first Brilloum zone (fc = 1), In this 
case, the system changes directly from the one-chain to 
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FIG. 10: The phonon spectrum for the two-chain configu- 
ration for different values of the density and with k = 1 
and a — 2. The transverse acoustical (TA), the transverse 
optical (TO), the longitudinal optical (LO), and the longi- 
tudinal acoustical (LA) phonon dispersion curves are identi- 
fied, respectively, by the solid, dashed, dotted, and dot-dashed 
curves. 



the three-chain regime. 

The dependence of the phonon dispersion curves on the 
linear density for the two-chain configuration, in the case 
with K — 1 and a ~ 2, is presented in Fig. 1101 In general, 
the frequencies of the longitudinal (transverse) branches 
increases (decreases) with increasing density. The expla- 
nation for such a behavior is similar to the one presented 
previously for the one-chain regime. 

The sound velocity, i. e. Vg — dio/dk |fe=o, corre- 
sponding to the LA phonon mode is shown in Fig. [TT] 
as a function of the density for a = 2. Notice that the 
sound velocity increases with density as expected. At 
n = 0.74155 a second order structural phase transition 
occurs to a two-chain configuration. At this phase tran- 
sition the monotonic behavior of Vg is changed because 
of the decrease in the density of particle along each chain 
and the continuous increasing of the separation between 
the chains. 

The general qualitative behavior of the phonon disper- 
sion curves presented so far for k = 1 is also found for 
different values of the screening parameter. 
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FIG. 11: The sound velocity obtained from the LA phonon 
mode for the one and two-chain configurations as a function 
of density for k = 1 and a = 2. 



V. CONCLUSIONS 

We studied a guasi-one-dimensional binary system of 
charged particles confined in a parabolic trap. The gen- 
eral structure and the main features of the normal mode 
spectrum were analyzed. The structural and dynami- 
cal properties were studied as a function of the density, 
the range of the interacting potential, which is an ex- 
perimental tunable parameter in systems like colloidal 
dispersions, dusty plasmas or even in a binary system of 
hard-spheres, and the parameter characterizing the bi- 
nary system, namely, the ratio between charges of the 
two types of particles. We found a very rich variety of 
ground state configurations, some of them are even not 
symmetric around the symmetry axis (i.e. i/ = 0) of the 
ID confinement potential. 

The number of chains as well as the internal structure 
in such chains are a function of the parameters of the 
system (a, n, and k). The set of structures was summa- 
rized in zero temperature phase diagrams which relate 
the ratio between the distinct types of particles and the 
screening parameter of the interaction potential with the 
linear density of the system. The structural transitions 
between the distinct phases were characterized as being 
of first or second order. A spontaneous symmetry break- 
ing in the charge distribution of the system was found 
corresponding to a structural transition that is charac- 
terized by a continuous change in the lateral position of 
the particles. A disordered phase was found in a small re- 
gion of the four-chains (case 2) part of the phase diagram 
where the particles are guasi— randomly distributed over 
three chains. The disordered three-chain configuration 
transits in a discontinuous way to the four-chain (case 2) 
ordered phase. 

The phonon dispersion curves for the one- and two- 
chain structures were obtained and analyzed for different 
values of a and n. In general, the normal mode frequen- 
cies depend on the linear density of the system. The 
transverse optical mode decreases monotonically with in- 
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creasing density, and such a feature was observed for all 
values of k. and a. The opposite behavior was observed 
for the longitudinal modes (optical and acoustical), i.e. 
there is a hardening of such modes with increasing den- 
sity. 

A remarkable softening of one of the branches of the 
dispersion curves was found at the second order struc- 
tural transition. The softening of the TO phonon disper- 
sion curve may occur for at the center or at the edge of the 
first Brillouin zone (k = it /2a), depending, respectively, 
if the transition is between the one- and the two-chain 
configuration or between the one- and the three-chain 



structure. 
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Appendix A 

The matrix w^I — D (where I is the unit matrix and D is the dynamical matrix) is used in the calculation of the 
normal modes for the one- and two-chains configurations. The dynamical matrix is of the form 

uj"^ - Ai -A 

Aa;2_^2 -Aq 

-As W^-Ag 

-^6 Aa;2 - A^ 

where Ao;^ = — Wq. The quantities Ax, A2 {A3, A4) are associated with the interaction between particles with 
charge qa (qb), while A^ and Aq account for the interaction between distinct charges. For the one-chain configuration, 
such parameters are given by 



Ai = a ^ 5 



3=1 



,-K(2j-l)/n 
(2j ~ 1)3 



2k(2j-1) , K2(2j-1) 



+ 



-2k j / 71 



4Kj , (2Kj)2 



2 + — + 



[1 — cos{2kja)], 



(Al) 



A2 = -aJ2 



-K{2j-l)/n 



1 + 



«(2j - 1) 



-E 



^ p-2Kj/n 



^ (2i)= 



1 + 



2Kj 



[1 — cos{2kja)], 



(A2) 



^3 = aJ2 



-K(2j-l)/n 



t (2J-1)^ 



2K(2j-l) , K2(2j-1) 



-2K,j / n 



t (2j? 



2 + — + 



n 



+ 



[1 — cos{2kja)], 



(A3) 



A4 = — a^^i 



-K(2j-l)/n 



U (2J-1)3 
-o?yn^ l-h 



/<2j-jO 



2Kj 



[1 — cos{2kja)], 



(A4) 
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^5 = —ay^n^^ 



-«(2j-l)/r, 



t (2i-i)^ 



3{k{2j - l)a)], 



(A5) 



^ p-K(2j-l)/n 



K(2j - 1) 



[cosik{2j - l)a)]. 



(A6) 



The dimensionless wave number k is in units of 7r/2a. 

For the two-chain regime the parameters in the dynamical matrix are of the form: 



'Kr/n 



^ — ^ r-^ I Kr"^ nr n"^ J n 



E 



4Kj (2«;j) 



2n 



[1 — cos(2fcja)], 



(A7) 



A2 ~ ay n 



— + — + 1 + — 

nr n 



2kj 



[1 — cos[2kja)\, 



(A8) 



^3 = «E 



3 3k 



(2j-l)^(- + - + -)-(l 



t (2j) 



n 



[1 — cos(2fcja)], 



(A9) 



A4 — a i 



r3c2 ^2^2 3^ 



KC , Kr , 



2 \ ^3 ' 



t (2jT 



2«:j 



[1 — cos{2kja) 



(AlO) 



i2j - 1)^(4 + ^ + 4) - (1 + -)! [--(fc(2j - IM 



(All) 



^6 = -"E ' 



-Kr /n 



— + ^ + (1 + —) [cos(^(2j - !)«)], 

r n nr n 



(A12) 



where r = -y/ (2j — 1)^ + c^, and c — d/a (d and a are indicated in Fig. [T]) is chosen in order that the eigenvalues of 
the dynamical matrix be positive. 
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